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We construct a coarse-grained effective two-dimensional (2d) hydrodynamic theory as a theoretical 
model for a coupled system of a fluid membrane and a thin layer of a polar active fluid in its ordered 
state that is anchored to the membrane. We show that such a system is prone to generic instabilities 
through the interplay of nonequilibrium drive, polar order and membrane fluctuation. We use our 
model equations to calculate diffusion coefficients of an inclusion in the membrane and show that 
their values depend strongly on the system size, in contrast to their equilibrium values. Our work 
extends the work of S. Sankararaman and S. Ramaswamy [Phys. Rev. Lett., 102, 118107 (2009)] to a 
coupled system of a fluid membrane and an ordered active fluid layer. Our model is broadly inspired 
' by and should be useful as a starting point for theoretical descriptions of the coupled dynamics of 

X^/^ \ a cell membrane and a cortical actin layer anchored to it. 
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I. INTRODUCTION 



Theoretical studies of equilibrium properties (both static and time-dependent) of fluid membranes are extensive and 
well-established by now [JQ. Mechanical properties of membranes are essentially controlled by a few parameters, such 
as the membrane tension, bending modulus and spontaneous curvature [l[ . These parameters completely characterise 
membrane fluctuation spectra at thermal equilibrium. By contrast, theoretical studies of the dynamical properties 
of driven fluid membranes are relatively at early stages of development. In this article, we consider the coupled 
I dynamics of a single-component fluid membrane and an active fluid layer anchored to it, within a coarse-grained 
hydrodynamic approach. Apart from its direct theoretical interests, there are biological implications as well. Studies 
of biomembranes are important in cell biology context, since all cells are covered by cell membranes which typically 
have complicated internal structures [|| . At a coarse-grained level, however, cell membranes are usually described as 
a bilayer fluid (with complicated microscopic internal structures, see, e.g., 0, H|). In eukaryotic cells, there is a thin 
layer of cross-linked actin filaments anchored to the membrane and associated with the cell cytoskeleton in typical 
eukaryotic cells. Knowledge of the properties of fluid membranes in equilibrium are not sufficient to characterise real 
biological membranes, such as plasma membranes of any cell, which are inherently nonequilibrium active systems 
0; 111,13] • For instance, cell membranes are constantly maintained out of equilibrium due to nonequilibrium processes 
^ of various types. It has now been realized that such nonequilibrium behaviour may underlie aspects of biomembrane 
dynamics [8|,l9(, previously attributed purely to equilibrium thermal fluctuations. These nonequilibrium aspects affect 
experimental measurements of physical quantities characterising cell membrane dynamics, e.g., one finds anomalous 
diffusion @, E3] (however, these are primarily due to coralling and inhibiting effects, different from active fluid effects 
\ discussed here). It is now generally accepted that such complex dynamical behaviour have their origins in generic 
7-H ■ nonequilibrium drives. In a living cell, there may be several different sources for such nonequilibrium behaviour - it 
could be due to active proteins inside the membrane that are often Adenosine Triphosphate (ATP) consuming enzymes, 
• i-h ■ or the membrane could be driven out of equilibrium by the nonequilibrium fluctuations of an active cortical actin layer 
anchored to the membrane. The former case has been studied both experimentally and theoretically in details, see, 
e.g., Ref. [ll|. The latter case is discussed in Refs. [Hl[II[. In particular, Ref. introduced a model that couples a 
membrane with actin polymerisation and contractile forces due to molecular motors, e.g., myosin. For a configuration 
of actin filaments which are typically normally directed to the membrane, Ref. [12[ has shown how motor activities 
lead to generic transverse membrane waves. Ref. fl3j discusses the influence of actin filament elasticity and retrograde 
flow on the force-velocity relations of motile cells. A related problem is the diffusion of small inclusions embedded 
in membranes, which is challenging both theoretically and experimentally. This has been extensively studied. For 
instance, Ref. [3| theoretically shows how curvature fluctuations (both thermal and non-thermal) in a membrane 
leads to enhancement of the effective thickness, that in turn causes a reduction in the measured diffusion coefficient. 
In a related work, Ref. discusses how membrane fluctuations, both static or quenched and annealed or dynamic, 
affect the diffusion coefficient. By using a mixture of analytical and numerical methods they find that for a membrane 
with quenched fluctuations, the diffusion coefficient is substantially reduced by the quenched roughness, and is nicely 
described by an area-scaling law proposed there. In contrast, for the annealed case with small bending rigidity, the 
reduction in the the diffusion coefficient is less than in the quenched case. These studies on diffusion coefficients, 
however, do not include any nonequilibrium effects that arise due to the active stress (see below) in the active fluid 
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layer anchored to the membrane. Elucidation of the generic effects of the active stress present in the ordered active 
fluid layer on the statistical properties of the attached fluid membrane is a major purpose of the present study. 

Although a cell cortex is typically not in an ordered state for the polarisation degrees of freedom (the actin filaments) , 
we consider an ordered state and study fluctuations about it. In particular, we assume the ordering direction to 
be parallel to the membrane. The other possible ordered state, viz where the actin filaments are predominantly 
perpendicular to the membrane is important as well (see, e.g., Ref. jll} for a study with perpendicular configuration 
of the filaments). Since a prime motivation of this work is to theoretically study in a general set up the effects of 
the nonequilibrium dynamics and fluctuations of an ordered layer of polar active fluid on a fluid membrane that 
is in contact with the active fluid layer, we have chosen to work with the parallel orientation as an illustrative 
example. Nevertheless, our scheme of calculations may be easily extended to filaments with perpendicular ordering. 
We briefly discuss this later. Our work here will have broad relevance and act as a starting point for more detailed 
study of the dynamics of a single-component fluid membrane in contact with a thin layer of cortical actin filaments. 
Recently, artificial actin layers attached to a fluid membrane have been constructed and experiments performed on 
them [la . [13]. Experiments include investigation of the structure of the actin layer and measurement of diffusion 
of inclusions in the membrane. The latter experiments do find rising lateral diffusion coefficients for larger systems 
[l6j | . While preparation of and experimental studies on artificial actin layers of a particular macroscopic configuration 
(i.e., a particular ordered state) anchored to a membrane is definitely technically a difficult task, some of our results 
described here may in principle be tested by performing controlled experiments on such systems. 

In this article, a thin cortical layer of actin filaments is considered at length scales much larger than the filament 
lengths and layer thickness with polar ordering, for which a generic coarse-grained continuum two-dimensional (2d) 
description would be appropriate. Further, we consider time scales larger than the unbinding time scale of the 
cross-linking proteins of the actin filaments, so that the actin network behaves like a fluid. In order to model the 
nonequilibrium dynamics of cortical actin layer, we use the active fluid description, proposed and elucidated in 
Refs. [H, [lj|. Such approaches have been successfully used in a variety of driven systems, see, e.g., Refs. [20l - [23j 
for extensive discussions about the subject. There have been a number of studies on the properties of active fluid in 
confining thin geometries before. For example, Refs. 0, HH proposed a mechanism for development of contractile 
rings in cell cortex. Ref. [25[ discussed large scale flows of the acto- myosin cortex in Caenorhabditis elegans zygotes. 
Further, Ref. [2(| provides a one-dimensional (Id) model for pattern formation in active fluids (see also Ref. 27] for a 
Id model of patterns in a cell cortex). In addition, Ref. J28| points out the existence of universal dynamic patterns by 
lateral membrane waves in motile cells. Our results here are complementary to those mentioned above. In this article, 
we systematically derive a set of effective 2d coarse-grained hydrodynamic equations of motion for the membrane 
height field, local orientational order parameter, concentration and the velocity fields, which we use to obtain results 
characterising fluctuations and correlations in the system. In particular, we show that there are generic moving or 
static instabilities (i.e., with or without propagating modes). We use a stochastically driven versions of these to 
illustrate the effects of the active (nonequilibrium) dynamics on the lateral and rotational diffusion coefficients of an 
inclusion in the membrane. We find that under certain situations there could be non-trivial system size dependences 
of the diffusion coefficients in comparison with their equilibrium expressions. The rest of this article is organised 
as follows: In Sec. |TT] we describe the basic model considered here and set up the equations of motion. We then 
show in Sees. IIII Al and 1111 Bl that these equations, upon linearisation show emergence of linear instabilities when 
the strength of the nonequlibrium drive, denoted by A/x exceeds a critical value. Finally, in Sec. IIVI we calculate 
lateral and rotational diffusion coefficients of an inclusion in the membrane; we show that for a membrane-active fluid 
combine lateral and rotational diffusion coefficients depend on the system size in ways that are very different from the 
known results for diffusion in systems at equilibrium. We summarise and discuss our results in the context of existing 
experimental results in Sec. |Vj 

II. MODEL SYSTEM AND EQUATIONS OF MOTION 

For the sake of simplicity, we analyse here the effect of activity on the fluctuation spectrum of a planar (flat on 
an average) membrane only. Our calculational framework as we formulate below is general. However, we use it to 
study, within a linearised approach, typically the properties of perturbations about an ordered, uniform reference 
state of the active fluid. The construction of our model is generally inspired by the structure of a real eukaryotic 
cell. There is a thin layer of cortical actin or cell cortex anchored to the cell membrane in an eukaryotic cell. This 
layer is rich in actin filaments which may be locally preferentially orientated either parallel or perpendicular to the 
membrane. The bulk of the cell is the cytoplasm consisting of cytoskeleton and other cellular organelles. We model 
the cortical actin layer by a thin film of active fluid of viscosity rj, thin in the z-direction and spread along the xy 
plane with active polar particles (since the actin filaments are polar) suspended in it. In our discussions below we 
only consider a macroscopically ordered state of the active particles (actin filaments) with in-plane ordering. We do 
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not discuss the other possible case of perpendicular orientations in details here. The active fluid is assumed to be 
anchored to one side of a fluid membrane characterised by a bending stiffness. We consider two distinct versions 
of our system: (i) Model I: The fluid membrane-active fluid combine is embedded inside a bulk isotropic fluid on 
both sides. The bulk fluid viscosity rj' is assumed to be much smaller than the active fluid viscosity (for simplicity; 
see below), and (ii) Model II: The system rests on a solid substrate below. The two cases are physically different: 
First of all, the presence of the solid substrate introduces friction and, as a result, the momentum (same as the 
hydrodynamic velocity v for an incompressible system) is no longer a conserved field. In contrast, when there is a 
bulk fluid surrounding the system there is no friction at the interface and hence v remains a conserved variable. This 
leads to long-ranged hydrodynamic interactions in such a system. Secondly, the presence of the solid substrate breaks 
the full three-dimensional (3d) rotational invariance of the problem, where as, when there is an embedding bulk fluid, 
the system is invariant under the full 3d rotational invariance (see below). This has important ramifications on the 
possible structure of the membrane free energy, as we will discuss below. Despite the cell biological inspiration for 
our work, the connections between our model and the structure of an eukaryotic cell is not very strict, and thus our 
model and results cannot be directly applied to a biologically relevant system: First of all, the cortical actin layer 
does not exist in an ordered state with in-plane orientation; at the scale of the cell it is generally isotropic. Secondly, 
in a real cell, there is not any sharp dividing surface separating the cortex from the bulk. Although our idealised 
model is not accurate for a real cell, nevertheless it allow us to obtain interesting and non-trivial result in a simple 
set up, highlighting how the environment (in the form of a solid substrate below or a passive, isotropic embedding 
fluid around the system) may drastically affect macroscopic physical properties. Our work may be considered as a 
starting point of more realistic calculations. Our treatment generalises Ref. [2!| to the case of a coupled system of a 
fluid membrane and an active fluid film. (See also Ref. [3(| for discussions on the effects of a solid substrate on the 
dynamics of thin liquid film covered by a membrane.) We now discuss the constructions of Model I and Model II in 
some details below. 



Consider first our Model I, where the fluid membrane-active fluid layer combine is embedded by an isotropic fluid 
on both sides. In order to keep the ensuing algebra tractable and to obtain closed equations and non-trivial results 
without having to introduce too many details we make a number of simplifying assumptions which we now discuss. 
First of all, we assume r\ S> rf. To set up the context for our work, we refer to the well-known results in Ref. [3lj ]. 
where the translational diffusivity of an inclusion is calculated in a thin 2d inflexible flat layer of an isotropic viscous 
fluid of viscosity n confined between bulk isotropic viscous fluids of viscosity r/i (rji/v finite) on both sides. We here 
consider a thin active fluid layer containing orientational degrees of freedom in their ordered state and covered by a 
fluid membrane on one side of finite stiffness. However, we consider only the special limit of n 3> 7]'. Thus our work 
here can be thought of as generalisation of and complementary to Ref. [31| • 



In order to specify the problem completely, we impose the following boundary conditions: When the fluid membrane- 
active fluid layer is covered by an isotropic fluid on both sides, the boundary conditions are as follows: (i) At the 
interfaces (z = hi and h 2 ) we impose p to be parallel to the local tangent plane on the surfaces, i.e., p • N = at 
z = hi and h 2 , as the current of active particles is along p and the particles cannot leave the film, where N is the 
local normal at hi and h 2 , and (ii) continuity of the shear stress at z = h\ and h 2 . The free energy functional T p of 
the system is a functional of hi, h 2 and p. The form of T v may be inferred from symmetry considerations. It must 
generally be invariant under an arbitrary tilt (equivalently a rotation) hi t 2 — > hi l2 + a • x of the free surfaces, where a 
is an arbitrary 3d vector and x is a 3d radius vector. This ensures that the most leading order (in gradients) coupling 
bilinear in p and hi or h 2 could be of the form V • pV 2 /ii.2- This is a polar term, since it has no p — >• p symmetry. 
Further, it violates h\_ 2 — > —h\_ 2 symmetry as well, which is admissible since the actin filaments are anchored only 
on one side of the membrane. For the polar order parameter, we use the Frank free energy [32[ in the limit of equal 
Frank's constants, denoted by D here. Assuming that the surface tension of the membrane is negligible, the general 
form of the free energy functional of the combined system of membrane (top) , free surface (bottom) and active fluid 



A. Model I: System covered by isotropic fluid on both sides 
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dz C(V-p)(V : 



■ 2 hi)S(z - hi) + C(V • p)V 2 M(z - h 2 ) + D{V lP3 ) 2 , 



(1) 



where k is the bending rigidity of the membrane, a the interfacial tension of the surface at h 2 , C is a coupling constant 
that couples the orientational field p with the fluctuations of the height fields. In principle, one needs to solve for 
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the dynamics of hi and h% separately for a full dynamical description of the problem. Operators Vi = {-§^t-§^) and 
V = (Jj, ^j, Jj) are the 2d and 3d Laplacians, respectively. At this stage, in order to simplify the ensuing algebraic 
manipulations, but still be able to display the dramatic effects of the hydrodynamic interactions, we consider a special 
simplifying (but admittedly artificial) limit of the problem with a being large enough (formally a diverging a), such 
that fluctuations of /12 are suppressed and may henceforth be ignored. Thus we are required to solve for the dynamics 
of hi only (along with p and c). It must be mentioned that the limit of large a is primarily of theoretical interests 
and is not expected to be observed in any real (biological) system. Despite these limitations of our assumption of 
large a, we obtain interesting results which are expected to be present qualitatively even with finite er. In order to set 
notations simpler, we set hi = h and /12 = in our subsequent analysis below. Thus hi drops out of the dynamics; 
see Fig. ([T]) for a schematic picture of our Model I system. Consequently, the free energy (TTJ) reduces to 
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FIG. 1: A schematic diagram of our model system showing the membrane and the active fluid layer. The arrows indicate the 
direction of macroscopic orientation (here along the :r-axis). 



?L[h,p] = \J d 2 r [n{V\hf] 

+ i J ' d 2 r J dz C7(V ■p)(V 2 h)5(z-h) + D(V lPj ) 2 . (2) 

It is invariant under h — > h+& ■ x. Now imposing the kinematic boundary conditions [33| on v z at h for an impermeable 
membrane we write: 

dh 

-^+v ± -V ± h = v z , (3) 

which connects the height fields h with the hydrodynamic velocity field v = (v x , v y ,v z ) at the location of the membrane 
(z = h). Here, vj_ = (v x ,v y ). In the expression fl2J) above, when D — C 2 / k > 0, the equilibrium phases are spatially 
uniform, else modulatedphases are possible. In this article we consider D — C 2 / k > 0, corresponding to uniform 
equilibrium phase only [34J . Thus finally the relevant slow modes in this effective 2d problem are (i) a concentration 
of active polar particles c(r, t), (ii) a local orientation field (describing the local orientation in the cortical actin) 
given by p(r, t), and (iii) a local height fields h(r,t) to describe the nearly flat membranes (in the Monge gauge [35j), 
where r = (x, y) is the in-plane coordinate and t is the time. We take the xy-plane as the easy plane for polarisation 
fluctuation. Further, we consider macroscopic orientational order given by a reference state p = (1,0,0) i.e. we have 
macroscopic ordering along the x axis and impose fixed length constraint on p: p = 1. For an active system, the 
polarity implies a current wocp with respect to the fluid, where vo is a characteristic drift velocity. 
The relevant intrinsic stress field of the active particles is given by [2(| H2, 



4 = A M c(r)p i (r)p i (r), (4) 

which is of nonequilibrium origin. Stress afj is said to be contractile or extensile for the constant A/i < or A/i > 0. 
The generalised Stokes Equation for v, is obtained by using the force balance condition after neglecting inertia as 
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appropriate for small masses in typical biological systems. Velocity components vj_ are to be solved from the full 
Stokes equation 



where a°^- 



N 



( a xji a yj)- At the membrane (z = h) we have p 

using the Monge gauge for the membrane we have N = (-VxM)/\/l 
membrane surface. This gives us 



0, (5) 
where N is the local normal at h(x,y). By 



Pz 



d x (h) 



(V±h) 2 to be the outward normal to the 



(6) 



This then yields d z p z ~ j-d x h. In order to solve for v from the above equations we must first find out V ■ a a . Within 
our linearised treatment, assume = x + 9y with small where p± = (p x ,p y ) ~ (1,9) to the lowest order. With 
\6\ -C 1 considering small fluctuation, we may obtain different components of Vj<7?-. To eliminate II we solve the 
Stokes equation for v z in the lubrication approximation, i.e., we write d z Tl = —dia? z = —A^i^-d 2 h. This yields 



U(x,y,z,t) = P Q + A^(h 



z )d 2 x h-f(h), 



(7) 



where f(h) = —^jjf- — —nV^h + C*V^V^ • px is the elastic force of the membrane. Equation ([7]) tells us that the 
balance of the normal component of stress at the membrane gives the condition of balancing the fluid stress by the 
clastic force of the membrane © . Thus the active contribution to the pressure in ([7]) comes with the same signature 
of A/i as in the active stress expression (QJ. 

The solutions of equation ([5]) are facilitated greatly in terms of the in-plane Fourier transformation: For a system 
as above, the general form for the hydrodynamic kernel is l/(r]q 2 + rj'q/ho) where ho is the average thickness of the 
system, where rj and 7/ are the viscosities of the active fluid and the surrounding passive isotropic fluid, respectively. 
Thus for q > r)'/(r]ho), the hydrodynamic kernel is identical to that of a free standing system. If (for the sake of 
simplicity) we now assume that r\ 3> rf, then in a broad window of wavevector, the hydrodynamic kernel may be 
approximated by 1/r/q 2 . Further in that limit, the continuity of the shear stress at z = h and z — actually implies 
vanishing of the shear stress at z = h and z = 0. Hence, in order to solve for v x and v y , we now need to solve 
Eq. (0 subject to boundary conditions of zero shear stress at z = h and z = 0. To the leading order in smallness this 
translates into 



at z = h and 



az oz 



r]— — = 



dz 



(8) 



(9) 



at z = (since p z = strictly at z = 0). Thus, for this wavevector range, effectively we have a free standing active 
fluid, covered on one side by a fluid membrane and a free surface on the other. In what follows below we stay in this 
limit only j36|. Thus 
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(10) 
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(ii) 



for qho <C 1 (we work in this limit which enables us to continue using the Lubrication approximation 
performed an in-plane Fourier transform with q as the Fourier conjugate of r. Functions ff,i 

(L{„ r ,\fL{ „ , ,\\ _ 2K B T r 



Here we have 
x, y are zero- mean, 



Gaussian noises with variances (//'(q, u))fj(— q, — u))) = ^ Sij. We have considered thermal noises for simplicity, 
although in real biological situations there are non-equilibrium noises as well. 

From the incompressibility of the fluid we obtain v z = — J Vj^ ■ ~v±dz together with the condition v z = at z = 0. 
Neglecting terms with higher derivatives and linearising about (h) = ho and (c) = cq and denoting h and c as the 
height and concentration fluctuations from their respective averages we obtain by using Eqs. (|10j) and (jlip 



dh 

m =Vz{z 



1 3Afj,c hl 2 2 

h) = — i — q qJl 



iChoq' k q y 9 — nhoq 6 h — (A/ihoq x c 



2c h A^,q x q y 9 + 2A^c ^/i)] + 6i 



(12) 
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in the Fourier space for a free standing system, £/j is a zero-mean Gaussian white noise related to ff 1 and hence with 

1K B Thl 

a variance — - — - . 

n 

Dynamics of the polar orientation field p differ from that of the more usual nematic director in that the equation 
of motion of p now must include terms [371 ] which violates the p > — p symmetry of the nematic director field. 
Apart from the usual terms [371 ] we can have a symmetry-permitted spontaneous splay term [29] in the free energy 
functional like J- sp = — J d 3 xCV • p, where C is a parameter which depends on the concentration of activity i.e. 

C = C(cq) + C'(cq)Sc + = C + C'Sc. Such terms may arise essentially due to the generic structural differences 

between the head and tail of the polar molecules. These will contribute — = — TC'Vc to the equation of motion 
for p, r being a kinetic coefficient. The 3d equation of motion for p y = 6 becomes 



d t 6 = -cuvodj - £d y c + CV 2 d y hS{z - h) + (XA yx - Q yx ) + DV 2 9 + & 



(13) 



where £ = TC', Ai 



2 ( V ^i 



is the strain rate tensor and fly = \{Si v j — VjVi) is the vorticity tensor. The 



first term of Eq. (fT3|) represents the advection of polar particles with the embedding fluid. The coefficient a\ is not 
necessarily unity, due to the lack of Galilean invariance. The next two terms are coupling terms coupling 6 with the 
gradient of concentration of the active particles and the fluctuation of the membrane surface respectively. Coupling 
constant A couples flow (strain rate tensor) with the local orientation. For stable flow- alignment |A| > 1 [32]. The 
last term is the diffusion term representing the Frank free energy (2(| contribution to the dynamics in the limit of 
equal Frank's constants. Noise £g is zero- mean, Gaussian distributed with a variance 2KbT [38j | . After eliminating v x 
and v y by using Eqs (jlip and (|10|) . under z-averaging yields an equation for 6, given by equation (|35p . as given in the 
Appendix (see Sec. IVI[) which depends only on h, 9 and c. When there is a solid substrate below, 8 follows the same 
equation as (j 13[) . except the tilt-polarisation coupling term now reads Cd y h. 

In our system the drift velocity of the polar particles with respect to the fluid is vo and hence the current due to 
active particles is given by vocp. Thus the continuity equation apart from diffusion is given by dtc+ V- ](v+uop)c] = 0. 
Now linearising about c — cq + c and using incompressibility we get 



d t c = -iv coq y 6 ~ ivoq x C + 0(q^ 



2\ 
,Qy)- 



(14) 



B. Model II: Membrane-active fluid combine rests on a solid substrate 



active fluid 



membrane 



h(x,t) 




solid substrate 

FIG. 2: A schematic diagram of our model system showing the membrane and the active fluid layer. The arrows indicate the 
direction of macroscopic orientation (here along the a;-axis). 



We now briefly discuss setting up the equations of motion when the fluid membrane-active fluid layer combine rests 
on a fully flat solid substrate; see Fig. ([2]) for a pictorial representation. Model II is important in the context in cell 
locomotion on a rigid substrate and lamellipodium. A lamellipodium is a cytoskeletal projection on the mobile edge of 
a cell. It is effectively a quasi-2rf network of actin filaments and is responsible for cell movement along a solid substrate 
(see, e.g., Ref. [|| for a general description of this). We essentially follow the line of arguments and calculations of the 
previous section, differing only in some details. First of all, in the present case we are concerned with the dynamics 
of the top fluid membrane (at z = h) - active fluid layer combined, where as the bottom solid surface (at z = 0) is 
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fixed. The latter is the preferred frame of reference for this model (hereafter Model II). This breaks the 3c? rotational 
invariance discussed for Model I above. This consideration then dictates the form of the free energy F$ as 



Fs[h} = ~ J d 2 rn{\/\hf + ~ J d 2 r J dz Cp ■ VhS(z - h) + D(V lPj ) 2 , (15) 

where, as before in n is the bending rigidity of the fluid membrane. Note that the bilinear coupling term in h and 
p is more relevant (in a renormalisation group sense) than the corresponding terms in ([2]) . We continue to denote the 
coupling constant by C. Further we use the equal Frank's constant limit, represented by D, for the Frank free energy. 
The pressure II may be simply evaluated as above; we obtain 

n = A/j l ^(h-z)d^h + KWih-CV± -p_l, (16) 

where we have used the boundary condition II = at z = (setting the zero of pressure) and balanced the normal 
component of stress at z = h. Velocity components v x and v y may be calculated from the Stokes Equation. For 
no-slip boundary conditions for v_l at z = 0, which we impose here, Stokes equation for vj_ may be written in the 
lubrication approximation (|V^v^| << |9 z v^|) 

v d 2 z v ± - v ± n - v j( x a Xj = o. (17) 

Now use (fTT)|) for II in Eq. ([17]) to calculate v x and v y : Equation (|T7|) may be twice integrated with respect to z 
to obtain v x and v y . The ensuing constants of integration are to be eliminated by using the boundary conditions 
Vx — = v y at z = and d z v y = 0, d z v x + A/zp z (rj_, z = h) = at z = h. The resulting equations of motion for 
v_l = {v x , v y ) read 

vi = A M g hz^j VAhdlh) A M g - h 2 z^j V x d 2 x h + ^ (y - fe) VxVi* - 

C ( z 2 \ Am ( z 2 \ 

' ' hz) V±d y 6 + — - hz) [{d x c + c Q d y 6 + coh^d^e^ + c o d x 0e y ] . (18) 



77 V 2 ) 7?V 2 
Equation of motion of h may be obtained by using the incompressibility of the velocity (just as before) . We obtain 

dh hi 5A/j/igc 2 2 3c A/j 2 4/t/i 6 .C/io 2 

at = "? [ ^ — fffl.^--4-fl.fc-— ff^ + »-3- ff*» 

-1/3(Am/i ^c + 2A^c h q x q y 9 + Afic q 2 x h)} (19) 

The equation motion for the polar orientation field is the same as for Model I above. After z-averaging one obtains 
the final explicit form given by Eq. (|36|) as given in Sec. I VII The concentration field obeys the same equation (fl4|) as 
for Model I. 

III. DYNAMICS AND INSTABILITIES 
A. Results from Model I 

Having set up the effective 2d coupled equations of motion for h(q, i), #(q, t) and c(q, t), we now examine the linear 
instabilities about a uniform state h = 0, = and c = 0. We begin with the extreme case with immotile but active 
polar particles, i.e., vo = but A/_t ^ 0. Clearly, in this limit, concentration c decouples from h and 9. We analyse 
the mode structure in the low q limit. Notable characters of the mode structures are 

• There are O(q ) anisotropic contributions proportional to A/zco/r?, which are of purely active origin. In detail: 
Assuming the fields to have time-dependence of the form exp(Ai), we have for the eigenvalues A of the linear 
stability matrix 



(A - l)A//cogg + (A + l)Afic g 2 x _ Aficpg 2 ± 1 (A - l)A/xc gg 
Ar]q 2 Ar/q 2 r\q 2 2 2r\q 2 



+ (A + 1)W _ }2 _ 4^4!^I { (A - IfeJ - (A + (20) 
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Eigenvalues are anisotropic (in the Fourier space) and have complicated dependences on q x and q y . It may 
however be noted that when 2q 2 — (A — 1)^/2 — (A + V)q 2 x j2 have a definite signature (either positive of 
negative) and (A — l)q 2 — (A + l)q 2 small, one may expand the square root in ([20|) binomially. One of the 
solutions of A is then proportional to A/i[(A - l)q 2 y - (A + l)g£]/[2q£ — (A — l)q y /2 - (A + l)q 2 /2] to the leading 
order in smallness, and is thus unstable just above or below the line (A — l)q^ = (A + l)q x , depending upon the 
signature of A/i. This is clearly indicative of instabilities for either sign of A/i. The unstable eigenvalue gets a 
particularly simpler form for large |A| ^> 1. We find one of the eigenvalues 

Afe,^)~^p^ 2 -a, (2i) 

displaying clearly the instability for either signature of A/u clearly. Further, for q x = and q y = one may 
separately obtain the eigenfrequencies A as 

A(q x = 0, q y ) = ^— - A/iC , 0, (22) 

A(q x ,q v =0) = -^,^±Ia M c . (23) 
V 277 

Notice that the mathematical origin for the instability for either signature of A/i lies essentially in the mixed 
or cross-coupling terms (i.e., the ^.-dependent term in the ^-equation and viceversa) and the fact that, for the 
leading order in wavevector q, the terms on the right hand side of Eq. (1121) have the same sign; similarly the terms 
on the right hand side of Eq. (|T3]) have the same sign. Let us try to see why this is so: There is only one source 
of activity in the model, viz, the active stress given by Eq. Now, this contributes to two different quantities 
in the effective 2d descriptions: (i) 2d active pressure, that depends on h and (ii) 2c? analogue of the bulk 3d 
active shear stress that depends upon 9. The dynamic of h and 9 depend upon both of them. This explains 
the presence of the cross-coupling terms in the dynamics. Because of the general structure of the Stokes Eq. 
and our approximation of small v z (lubrication approximation), the 2d active pressure and the 2d shear stress 
come with the same signatures, the feature that is responsible for the occurrence of the instabilities mentioned 
here [39j . Since physical fields h(r) and 9(r) depends upon the Fourier modes h(q) and 0(q) for all q, h(r) 
and 9(r) display instability for both signs of A/i, at short enough wavenumbers. Note that these g-independcnt 
contributions do not mean the modes have finite life time at vanishing q; such g-independent behaviour is 
nothing but a consequence of using the Stokes' equation (instead of the Navier-Stokes equation) for the velocity 
field. In the Stokes' approximation one neglects the inertia in comparison with the viscous term. Of course, for 
very low wavenumber this is no longer valid and inertia effects will be important. Thus the above expressions of 
A cannot be used in the limit q — > 0. Note, however, that this does not question the mathematical validity or 
consistency of the Stokesian hydrodynamics that we have used. At the linear level, all modes labeled by q are 
decoupled. Hence, the properties of one particular mode is unaffected by any other mode with any wavevector 
(including vanishing wavevectors) . Since the initial amplitudes of the (small) perturbations considered here can 
be taken to zero smoothly, the validity of the Stokesian hydrodynamics with linearised approximation remains 
sound. 

At the 0(q 2 ) the contributions are easily presentable in the limits q x — and q y — 0, when their expressions 
simplify considerably, the contributions are still purely real and display instability just like the contributions at 
q°. Parts of the contributions are active, parts are however from the equilibrium part of the dynamics. We find 

A{q x = 0,q y ) = hzl^^-D^O, (24) 
M n\ 2A/iC . A/JCq/Iq 2 A+ 1 2 

Hqx,q y = 0) = — + — ^ — q x , -^-A/iCo - Dq x . (25) 

A schematic diagram of the eigenvalues is given in Fig. [3l We thus see that there are both stable (coming from 
D > 0) and unstable active contributions at 0(q 2 ). Thus depending upon the numerical value of A/i, there may 
be instabilities at 0(q 2 ). Ignoring anisotropy, Eq. (f24|) or (|25|) allows identify a length scale 

L ' ^ JM- (26) 



such that systems with linear dimensions L > L c are unstable. Note: Due to the anisotropy, the length scale 
L c should depend upon the polar angle which may be obtained in a more detailed calculation. Propagating 
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FIG. 3: Schematic plots showing the eigenvalues A in the two limits up to 0(q 2 ) 
for Model I for A/i > - left: A(q x = 0, q y ) vs q y (a downward parabola, continuous line) , right: A(q x , q y — 0) vs q x 
(upward (dot-dash) and downward (continuous line) parabolas). The broken horizontal line in both the plots denote A = 0; 

its meeting point with the downward parabola yields q c ~ 1/L C . 



modes due to the tilt-curvature coupling appear at even higher q, at 0(q 5 ^ 2 ). These propagating modes are 
strongly damped at smallest wavenumbers. Finally, at a sufficiently large q, all the instabilities are cut-off by 
the bending stiffness k of the membrane. 

What could be a typical value for L c l We take linear size of an active particle ~ lnm, volume fraction <f> = coa 3 = 1 
(since we assume to be in an ordered state), Dr) ~ K ~ 10~ e dyne/cm where if is a Frank elastic constant. Estimation 
of A/i is more ambiguous: We use the fact that approximately IKCal energy released per mole of ATP due to its 
hydrolysis. Since 1 molar mass of ATP ~ 500, we obtain from its definition A/x ~ 7KCal/(500gm/10 23 ) the free 
energy release per unit mass per molecule. All these, however, lead to L c ~ lOnm, a value rather low compared to 
the linear dimensions of an eukaryotic cell. However, we have ignored anisotropy while estimating L c . Moreover, for 
a realistic situation the effects of the ambient fluid is likely to be significant and should affect L c . In any case, a direct 
comparison with a living cell is not very effective due to the simplifying approximations and the assumption of an 
ordered state of the active particles that we have made. 

If we now consider motile particles (vq ^ 0) but continue to ignore c, the most dominant effect is that one of the 
modes will now pick up a propagating part with speed vq and dispersion proportional to q x . Although this is still 
subdominant to the most leading order instabilities, for sufficiently large vq, the propagating modes will be observed 
and the instabilities at higher q corresponding to this particular mode will now be moving. Finally, when the dynamics 
of c is considered, the eigenmodes have complex forms as a function of q. A simple way to analyse them is to look 
at special limits of either q x — or q y = (an alternative way, considered in Ref. [29] . is to examine the limit 
TC/vq > 1). We find 



-4« co^] 1/2 , (27) 

2Afic hl 2 nh 4 2A^c 

HQx,q y ^0) = q x q x , 

irj rj rj 

(A + 1) - Dq 2 x - iaxv q x , -iv q x . (28) 

Thus again in the intermediate wavenumber range there are 0(q 2 ) instabilities for sufficiently large Afj,. In addition, 
there are propagating waves coming from the concentration-orientation coupling. Thus the instabilities and the 
resulting patterns will be moving. Finally, at even larger wavenumber, all eigenmodes are stable owing to the bending 
stiffness of the membrane. If A/i = 0, the system is stable which is expected in the equilibrium limit. In Fig. 3] we 
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pictorially show the presence of instability for both signs of A/j, in plots of vs D and Afic Q vs k for eigenvalue 

^■{ixily = 0) for a given value of q x . The line in the D — A/ico/r] plane are given by Dq% = (A + l)Afxco/2rj for a 
given q x and the line in the k — A/iCo/77 plane is given by nhoq^. — —2Afico + A/icohQq^/2 again for a given q x . The 
plots clearly indicate that for any of the signatures of Afi some regions in either of the plots display instability. 




o 




I 




FIG. 4: A schematic diagram showing the unstable and stable regions in a phase space diagram for Model I with a given q x 
(left: A ^ CQ vs D , right: A/ico vs k) for eigenvalue A(q x ,q y = 0). The lines are given by the zeros of A(q x ,q y = 0) = for a 
given q x , ignoring v . 



B. Results from Model II 

When there is a solid substrate, the long wavelength dynamics is same as that in Ref. [2^|. We present them here for 
direct comparison with the free standing case. The two eigenmodes A of the linearised coupled dynamical equations 
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of h and 9 have the form (setting vq = 0) 



24^ q " x YhT qy 6^ r 



1 llA/xcpfo 2 , 2 A/xcpfeg 2 2 AAMCpfe| 2 a ^ A/xcp/ig ^ 2 4rC^ 2 2 

: 2 [{ — 12^- - -^-^ - Dq + ~^r q } - 8tC ^r qx% ~ ~^r^ q 



13A M coJ>feg 2 2 13(AjO a cgf»* 13(A^) 2 cg^ 2 2 , ^(A^) 2 ^ 2 
As in Ref. (29l |. to the lowest order in wavevector, there is a traveling instability with a growth rate proportional to 

1/2 

QyQx ; the mode displaying an instability travels in the +x direction (-x direction) if A/i > (A/i < 0). If we now 
include Vq, but still ignore the dynamics of c, the eigenfrequencies in the limits q x = and q y — are 



12?y 16?; 3^ 

Gr/ 



.io 1 « bfe -D^ + (2A+l) A ^M^, (30) 



Afe=0, % ) = -^ + (2A-l)^^±i 



M-(2A-1)^*}»-^ 



(31) 



It is worth pointing out the basic differences between the eigenvalues for a free standing system given by (I24I25[) and 
those for a system with a solid substrate below given by ([29)) or those in Ref. (2{| , whose calculational framework has 
been largely used here. In the former case, eigenvalues are (/-independent for small q. This is a consequence of the 
long-ranged hydrodynamic interactions. In contrast, when there is a solid substrate, the eigenvalues smoothly tend to 
zero as q — ► 0. This is because any long-ranged hydrodynamic interactions are screened by the solid substrate below. 
This means, no matter how small the activity (A//) is, for a large enough system (with linear dimension L > L c ) the 
system in Model I gets unstable. In contrast, for model II the screening of the hydrodynamic interactions yields that 

as long as the system has a thickness smaller than the critical thickness h 0c ~ \J J^^^ for a given A/i, the system 

with arbitrarily large in-plane linear dimensions remains stable. This crucial difference apart, eigenvalues (130131)) 
show instabilities for sufficiently high A/j. Although the latter features are similar to that for a free standing system 
(Model I), the underlying physical mechanisms are different: Unlike Model I, where hydrodynamic interactions are 
responsible for type of the instabilities observed there, for Model II, the tilt-concentration coupling characterised by 
the coupling constant C [29| is responsible. 

Regardless of the details of the boundary conditions on the system, qualitative implications of both hold for a 
variety of phenomenologies. First of all, since v x and v y depend upon h, 9 and c, as given in Eqs. (|11[) and (|10[) for 
Model I, and Eq. (fTSJ) for Model II, as soon as one of the modes become unstable, the initial homogeneous non-flowing 
state will be unstable and the active fluid starts to flow. While our simple model calculations by themselves cannot 
capture the full phenomenologies of flows in cell cortex (see, e.g., Refs. 0, HH; see also Ref. [26| for a Id model 
for pattern formation in active fluids), our results open up possibilities of flows in the cell cortex within calculations 
in a simple setting. We show instabilities develop under generic conditions. In addition, there may be waves under 
various conditions. Our results are universal in the sense that they do not refer to any specific cell or do not depend 
upon very specific biochemical processes. Although our system and the results which follow cannot be directly related 
to any biological systems due to our several simplifying approximations, we believe our results will insp ire more 
realistic studies on (moving or static) patterns ultimately formed in cell membranes (see, e.g., Refs. [28l. l40l l4l|). 
Finally, in view of our results here, let us briefly consider the possible structure of the effective 2d model with a 
polar order that is perpendicular to the surfaces at z — and h. While quantitative predictions require detailed 
calculations, we can already make the following observations, based upon the framework developed above. In this 
case, again assuming a flat bottom surface (either a solid substrate or a free active fluid surface with a very large 
surface tension) the boundary conditions on the polarisation p arc p z = 1 at z = and p • V/i = 1 at z = h. For small 
fluctuations of p x and p y (fluctuations in p z will be second order in smallness) the boundary condition at z = h implies 
Px = 9 x h/2, p y = dyh/2. Further, p x and p y are constrained to be zero at z = 0. Thus, in a z-averaged description 
as above, p x and p v are slaved to Vft. Hence, staying within the framework of 2d effective description as developed 
above, the effective dynamics will be described by 2d equations of motion of height h and concentration c only. This 
effective dynamics is now expected to display signatures of instabilities and a spontaneous flow transition akin to the 
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Frederiks transition of equilibrium nematic liquid crystals at a critical thickness for contractile active stress with a 
given value of A/i, similar to those discussed in Refs. [H|,|4!|. Details of these will be discussed elsewhere. 



IV. DIFFUSION COEFFICIENTS 



The in-plane diffusive motion of membrane bound protein are a key ingredient in a great many biological functions 
including exchange of information, material and/or energy. A complete picture of how cells function and interact with 
their immediate surroundings requires understanding of this diffusive phenomena. The recent progress in experimental 
techniques to measure lateral (Dl ) and rotational (Dr) diffusion coefficients, like fluorescence correlation spectroscopy 
[42] ]. single particle tracking [43 |. or fluorescence recovery after photobleaching has revealed that many of the 
functions performed by proteins are crucially influenced by the diffusive behavior of the proteins (45j | . Apart from its 
obvious biological significance, diffusion of membrane-bound proteins is a good example of 2d diffusion in a fluctuation 
background (here the membrane), where the fluctuations are of nonequilibrium origin. Calculation of Dl in a strictly 
2d flow is subtle due to the divergence associated with the solutions of the 2d Stokes' equation, known as the Stokes' 
paradox in the literature pjjj . In order to overcome this Stokes' paradox, Saffman and Delbriick [3l| considered 
the mobility of a very thin, rigid object in a narrow almost 2d perfectly flat fluid layer that is surrounded on both 
sides by a further liquid and obtained finite results for Dl- Calculation of Dr does not suffer from any such subtle 
issues. Recently in Ref. [47| the authors calculated Dl for a protein molecule in a bare membrane. In none of these 
theoretical examples, the nonequilibrium nature of the membrane dynamics and the associated active cortical actin 
layer have been considered, although the cellular cytoskclcton is often linked with the mobility of a protein molecule 
in the membrane, see, e.g., Refs. [7J. We here elucidate possible effects of the cortical actin layer nonequilibrium 
fluctuations on the measured value of the diffusion coefficients within our effective 2d coarse-grained model. As we 
shall show below, the ensuing algebraic details is rather complicated. Hence we confine ourselves here to a treatment 
at the scaling level only. This suffices for our purposes here. 

Let us begin with the formal definition of the diffusion coefficients: Since we have an anisotropic system, we 
consider the general lateral diffusion tensor D^ which is defined by the relation (ri(t)rj(t)) = 2D^t, where r.i{t) is 
the Lagrangian coordinate of the protein molecule in the membrane at time t. Equating the Lagrangian velocity of 
the particle with the local in-plane 2d velocity field, away from the instability threshold the above definition leads 
to Df- = 5/ (27r) 2 ( v * fa' u ~ Qi 00 = 0)).We calculate the two lateral diffusion coefficients D xx and D yy below. 

Similar to the lateral diffusion coefficients, the rotational diffusion coefficient Dr is given by Dr — i J (\fl(q, u> = 
0)| 2 ), where fi is the z-component of the vorticity tensor defined above. Thus enumeration of diffusion coefficients 
requires calculations of certain velocity correlation functions which may be done by stochastically driving the 2d 
equations of motion developed above. In order to simplify our calculations we ignore any local effect of the protein 
molecule on n [47*]. Further, we set vq = 0, thus concentration c decouples from the dynamics at the lowest order 
in q. We present the main results below, with some of the algebraic details available in Appendix II. For a free 
standing system, we add zero-mean conserved Gaussian distributed thermal noises in Eqs. (110H11[) and a zero-mean 
thermal Gaussian distributed white noise in Eq. (f!J5j). Finally, we calculate the diffusion coefficients for immotile 
active particles, i.e., vq = 0. In addition, we ignore the concentration field c in the following calculations. 

We consider Model I first. Calculations are considerably simplified if we consider the dynamics only at the long 
wavelength limit, which suffices for our purposes here. In that limit, correlators of h, 9, v Xl v y are available in 
Appendix II. The noticeable feature is that all of Eqs. f4"Tj) have parts which diverge beyond a critical system size 
L c . Ignoring anisotropy, the length scale L c is given by Eq. (|26p above. From the correlators (|41l) in Appendix II. we 
obtain for the lateral diffusion coefficients (we show only the diverging parts below) 

KsT f gg C-^f 1 , ^(A-l) 



h J (27T) 2 [IV + (A-ffiA/. gg - ( A +y^ P 77 V 2q x 

KbT_ f jPq_ (^T) 2 1 / fr(A-l) V 

h J (2n)2 W + ^^ql V 1 V 2q x J 1 ^ 

Due to the complicated nature of the expressions Q32p we do not attempt to evaluate them exactly. Instead, we treat 
them at the scaling level. For a free standing system without activity, D xx and D yy depend on the system size L 
logarithmically, i.e., as logL, which is the equilibrium contribution. This contribution survives even when A/i =/= 
which we do not show explicitly above. Evidently, for an active system, the dependence of the active contributions 
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on L are very different, as given by Eq. (|26|) . Thus at the scaling level we obtain 

D xx ,D yv ~^Li(L 2 c -L 2 )- 2 . (33) 
n r] 

for L <C L c . Thus as L increases, D xx and D yy rise with L in power law fashion, unlike in equilibrium systems 
where such rises are logarithmic in L. Similarly, Dr, for a membrane on a thin non-active fluid film (i.e., in thermal 
equilibrium), has no dependence on L, it instead depends on the small-scale cut-off I (of the order of the particle 

size) as l~ 2 . In contrast, when the fluid film is active, Dr ~ jjj=f?yz below the threshold. Thus, similar to 

D xx and D yy , increases as L increases. For a thin inflexible isotropic fluid layer with viscosity rj surrounded by a 
bulk fluid of viscosity i]i with 771 /rj finite, theoretical calculations of Ref. [3l[ diffusion coefficient of an inclusion of a 
small but finite size is finite. Our work is thus an extension of Ref. [3l[, incorporating effects of orientation degrees 
of freedom and a fluid membrane, but staying at the special simple limit of 771 / 77 — > 0. In addition, our assumption of 
inflexible active fluid-isotropic bulk fluid interface will not strictly hold in a real cell. Thus more refined calculations 
are needed for better quantative estimations. We conclude this Section by mentioning that there are two important 
physical effects which we have not considered. First of all, it is well-known that for a membrane with a finite thickness 
in thermal equilibrium, fluctuations increase the effective thickness fl4| which in turn reduces the diffusivity. Similar 
effects should be observed in the present problem as well, whose quantitative enumeration requires further work which 
we do not discuss here. However, our results above holds in the limit of zero membrane thickness (the parameter 
ho here corresponds to the average thickness of the active fluid layer, and not of the thickness of the membrane at 
the top). Secondly, in view of the findings of Ref. [HI that a quenched rough surface reduces the effective diffusion 
coefficient substantially, whereas, for an annealed surface the reduction is relatively small, it would be im por tant to 
investigate the analogous effects in the present problem. However, direct application of the results of Ref. [15[ to our 
problem is difficult due to the presence of additional degrees of freedom (polarisation fluctuation 8) with long-ranged 
correlations and the active stress. Further work is necessary to settle this issue properly. This is beyond the scope of 
the present work. 

In order to calculate diffusion coefficients of an inclusion for Model II we need to set c = 0, since otherwise there 
are underdamped propagating modes in the system. From the correlators (|43l) we find that, unlike the free standing 
system, there are no critical lateral size beyond which instabilities set in; further their dependence on L is \ogL. 

Instead, now there is a critical thickness ho c such that for a film with thickness larger than ho Cl h$ c ~ y C ^A^ 



spontaneous flow instabilities akin to the Frederiks transition in equilibrium nematics set in |48l |49| . For a system 
with ho < hoc we find 



K B T h 4 
" r,ho (h 2 0c -hi) 2 



D XX ,D* ~ _Qc io g £ (34) 



for ho smaller than the critical thickness hoc- Thus D xx and D yy depend explicitly on ho, a measure of the system 
size. Note the differences between the expressions of diffusion coefficients in Model I and II. For Model I, they depend 
upon the linear size L of the system in a way markedly different from the log L dependence as observed in equilibrium 
systems. In contrast, in Model II, they depend on L as logi; in addition however they acquire non-trivial dependences 
on ho- 



V. SUMMARY AND OUTLOOK 



In this article we have set up 2d coarse-grained equations for a coupled system of a fluid membrane and a thin 
layer of active cortical actins anchored to it in terms of a height field, orientation field and concentration of active 
particles. We considered two cases of a free standing system and a solid substrate under the system, separately. The 
distinguishing feature of the former case is the presence of long ranged hydrodynamic interactions. We discuss the 
generic instabilities and patterns which appear due to activity. We use our equations to calculate lateral and rotational 
diffusion coefficients of an inclusion in the membrane. For the case of a free standing film, we have assumed that 
the surface tension of the bottom free surface is large, so that fluctuations of that surface are suppressed and are not 
considered in the subsequent calculations. This is mainly a theoretically interesting case and for real experimentally 
testable systems, surface tension of the surfaces should be finite and hence the bottom surface will have fluctuations. 
Despite the limitations of our approximations, our results bring out the differences between the two cases (free standing 
system and system resting on a solid surface) very clearly. For the sake of analytical convenience we have ignored 
a few details of cell membranes, e.g., effects of the surface tension of the membrane, the role of active proteins in 
the membrane and local modulation of the bending stiffness due to the inclusion; these may however be important 
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in a real biological set up. Qualitative features of our basic results should in principle be testable in standard cell 
biology experiments measuring, e.g., measurements of diffusion coefficients. However, direct comparisons with in vivo 
experiments will not be easy primarily due to our assumptions of idealised model and also due to many complicated 
features of a cell membrane (e.g., a cell membrane is actually a bilayer or the cortical actin is anchored to the 
membrane only at discrete junctions). Our work here relates to Ref. [3l| in that we include new effects coming due 
to orientational degrees of freedom, active stress and finite membrane stiffness, but consider the limiting case where 
the ambient fluid viscosity is much smaller than the active fluid viscosity. Despite the limitations of our simplifying 
assumptions, out work shows the dramatic effects of hydrodynamic interactions and provides a first step towards more 
realistic calculations. We hope our work will stimulate more realistic calculations in this direction. 

Apart from its phenomenological importance, we believe our work is a first step in formulating and understanding 
the 2d coupled dynamics of a fluid membrane and driven orientational broken symmetry modes. It would be interesting 
to see how the predictions of Ref. [50| get modified due to the fluctuations of the membrane. Secondly the fluctuating 
orientational degrees of freedom should create an effective long ranged interactions (of nonequilibrium origin) between 
different parts of the membrane. It would be interesting to see whether such interactions may allow a finite temperature 
crumpling transition of a 2d fluid membrane !51j, something which is prohibited in equilibrium (35[. 



VI. APPENDIX I: FULL EQUATIONS FOR 9 



In order to obtain an effective two-dimensional description of the dynamics of 9 for a free standing system (Model 
I), we average over the z direction to get 
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The corresponding dynamical equation for 9 when there is a solid substrate below (Model II) is 



d t 9 



-iaiVoq x 9 — i£q y c + iTCq y h — Dq 2 9 H — — - — -coAfihoq x qyh AfiCohQq v q^.h 



24?7 



_A_ 
16ry 



C 

+ '-^hlq x q y q i h - i^-h^q 2 q x 9 



Xk , 
V 

A ^„ 1,2 2a , X<7h 2 



(A-l) 
6rj 



A 



hgA^q x q y c+ - L Afic Q h 2 l q 2 9 - -^-c h^q y 9 



A/Lt 



,2 2, 



3r? 



Gr/ 



q q x q y h. 



(35) 
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VII. APPENDIX II: CORRELATION FUNCTIONS 



A. Model I: Free standing system 

Equations (fTTj) and (fTO)) reduce to 

i 

v x = —[c Q h Afiq y 9 + 2c Afiq x h], (37) 

v y = — \q x 9. (38) 

m 

The equations for 9 and h in the Fourier space in the long wavelength limit are 

iu>6 + Dq 2 9 - - — ^- (h c Afiq 2 9 + 2c Afj,q x q y h) - ^ + coA/xg^g = g$, (39) 
Zriq^ho y 2r\q z 

itoh \r(2c Q h Anq x q y 9 + 2Aficoqih) = gu, (40) 

m 

where gg are g^ are zero-mean Gaussian distributed white and conserved noises respectively. In the above we have 
ignored the dynamics of c and set v = 0. Equations (|4"0"j) and P0|) may be solved in a straight forward way to obtain 
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the corelators (\h(q,ui)\ 2 ) and (|#(q,u;)| 2 ). The resulting expressions are lengthy and not very illuminating; we do 
not present the full expressions here. Instead, we obtain the correlators in the long wavelength limit and use them to 
obtain correlations functions (\v x (q,uj)\ 2 ), (\v y (q_,ui)\ 2 ) and (|f2(q, ui)\ 2 ) (where fi(q) is the Fourier transform of the 
^-component of V x v). We find (ignoring parts which show no divergence) 



2\ _ Of_T/. LJM 
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Evidently, all the correlators diverge beyond a typical critical system size L c defined above. 

B. System in contact with a solid substrate 

The stochastically driven 2c? equations of motion of h, 9, v x and v y are 
= -i q x h-i — q y + £ x , 

. A/^Cq/Iq 



2>) 



dh UcoAfihlq 2 2h%coAnq x q y 

at ~ 12^ — h 3^ e + ^ 

dO 2 7(A - l)c h Afiq x q y \c Q h 2 Afig 2 A^c h 2 g 2 y Afic h 2 g 2 

m = ~ Dq9+ 24^ k+ 3^ ° 677 9 + ^ 

(42) 

Here the noises £ x , ^, and £| are all zero- mean Gaussian white noises. Equations (|4*2"j) lead to the velocity 
correlators 

.. . n . |3 . ^c^ ^ m^ - o)gj(-g," = o)) + {f|^m(q^ - 0)l 2 ) 

(|u x (q, a; = 0)| ) = ( — — v-*-) — ~- 

^ \Da 2 - XAflCoh o a 2 _ A Mcofeg 2 i coA/jfag (■ 7(A-1) , li o 1 

" 3?7 " 6t; " x I 39 6 J 



(K(q, W = 0)| 2 ) = ( 



-<£(q,w = 0)g(-q, W =0)), 

A M Co^g x ,2 <«fo> w - 0)C|(-q, W = 0)) + {^aL}a(|^( q , w = 0)| 2 ) 



2r/ 



n 2 _ AAj£Cojj 2 _ A^cpfeg 2 , coA^hg f 7(A-l) i-, 2 
^ 3)) y 6i) q x r) I 39 T 6 J y y 



+ (^(q,w)^(-q,w=0)> (43) 

Thus there are no instabilities at any finite lateral dimension of the system, and hence there is no critical L, unlike 
the free standing system. 
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